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ABSTRACT 


The  power  spectral  density  of  the  fading  power  of 
high-f r equency  radio  waves  and  methods  of  estimating  this 
function  are  the  subject  of  this  report. 

Evaluations  of  the  power  spectral  density  estimates 
obtained  by  two  different  methods  are  made.  The  asymptotic 
low- frequency  variation  of  estimates  computed  on  a  digital 
computer  and  the  expected  value  of  estimates  computed  with 
an  analog  computer  are  derived.  A  comparison  of  the  tv/o 
methods  is  made  and  power  spectral  density  estimates  com¬ 
puted  by  both  methods  are  shown. 

Power  spectral  density  estimates  are  shown  for  eight 
different  data  samples.  In  general,  the  power  density  is 
confined  to  frequencies  below  one  cycle  per  second. 
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1. 


Introduct ion 


The  purpose  of  this  investigation  was  to  make  a  study 
of  the  statistical  character  of  fading  radio  waves.  In 
particular  the  power  spectral  density  of  the  amplitude  of 
the  fading  waves  was  studied.  Several  samples  of  data  were 
statistically  analyzed  in  order  to  obtain  estimates  of  the 
true  power  density  spectrum  of  the  underlying  random  pro¬ 
cess.  The  estimates  were  computed  in  two  different  ways: 
the  first  way  utilized  a  digital  computer  and  the  second 
an  analog  computer. 

A  theoretical  evaluation  of  the  methods  of  computing 
the  estimates  of  the  power  density  spectrum  was  also  under¬ 
taken  in  order  to  aid  the  interpretation  of  the  estimates. 
The  expected  value  of  the  estimates  that  were  calculated 
with  the  aid  of  an  analog  computer  was  derived  and  the 
nature  of  the  variation  of  the  expected  value  of  the  digi¬ 
tally  computed  estimates  was  ascertained,  A  comparison  of 
the  digital  method  with  the  analog  method  was  made. 

The  next  tv/o  sections  of  the  report  deal  with  the 
nature  of  the  data  and  the  ionosphere.  The  fourth  section 
discusses  the  statistical  models  characterizing  the  observed 
process.  The  next  three  sections  describe  the  digital  com¬ 
puter  program  for  obtaining  estimates  of  the  power  spectral 
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density,  statistical  tests  of  these  estimates  and  the 
asymptotic  variation  of  the  expected  value  of  these 
estimates  for  low  frequencies.  The  last  part  of  Section  7 
contains  a  comparison  of  the  digital  computer  method 
and  the  analog  computer  method  of  estimating  the  spec¬ 
trum,  The  next  two  sections  describe  the  analog  com¬ 
puter  method  and  a  derivation  of  the  expected  value  of 
the  power  spectral  density  estimates  obtained  with  this 
method.  The  last  two  sections  present  the  results  and 
conclusions  of  the  study. 

2,  Nature  of  the  Records 

The  records  of  the  amplitude  of  the  fading  radio 
waves  were  obtained  by  recording  on  magnetic  tape  the 
detected  voltage  out  of  the  IF  strip  of  a  receiver.  The 
receiver  was  adjusted  so  that  the  output  was  directly 
proportional  to  the  radio- fr equency  input  over  the 
operating  range. 

For  some  of  the  data  the  signals  received  by  the 
receiver  had  been  transmitted  over  a  6000  km  path  and 
most  probably  were  reflected  three  or  four  times  by  the 
ionosphere.  Signal  frequencies  of  11.9,  35.1  and  17*9  Me 
were  used.  These  signals  were  c w,  keyed  for  i  lentifica- 
tion  every  half  of  an  hour.  The  remainder  of  the  data 
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consists  of  recordings  of  the  amplitude  of  WWV  at  20  me. 

Examination  of  data  samples  (Figures  1  and  2)  shows 
that  the  signal  underwent  considerable  variation  but  that 
the  variation  was  not  rapid,  usually  requiring  more  than 
one  fifth  of  a  second  for  any  considerable  change  in  mag¬ 
nitude. 

3 .  Description  of  the  Ionospheric  Phenomena 

The  mathematics  describing  the  statistics  of  the  am¬ 
plitude  of  the  reflected  radio  wave  should  be  based  upon 
some  physical  model.  The  physical  model  envisioned  in 
most  attempts  to  explain  the  fading  phenomenon  is  an 
ionized  medium  which  does  not  have  a  uniform  ionization 
density  throughout  its  extent.  This  medium  is  called 
the  ionosphere.  Due  to  the  non-uniform  ionization  some 
regions  of  the  ionosphere  will  reflect  radiation  more 
effectively  than  other  regions.  These  regions  of 
ionization  are  continually  forming  and  dissipating  and 
are  moving  with  different  velocities.  The  motion  of 
these  regions  causes  radio  waves  to  suffer  a  Doppler  fre¬ 
quency  shift.  The  radio  wave  that  is  finally  received 
from  this  scattering  process  then  is  a  v/ave  which  might 
be  thought  of  as  being  composed  of  components  differing 
slightly  in  frequency  and  differing  in  propagation  path, 
which  combine  in  such  a  way  as  to  yield  a  wave  varying  in 


amplitude  and  in  phase. 
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The  treatment  of  the  problem  depends  upon  several 
parameters.  These  are  (a)  the  size  of  the  scattering 
regions  (* spots’)  in  the  ionosphere,  (b)  the  number  per 
unit  volume,  (c)  the  rate  of  formation  and  disappearance 
of  these  regions,  (d)  the  average  velocity  of  these  spots 
and  (e)  the  uniformity  of  motion  of  these  spots. 

Whether  many  or  few  scattering  regions  in  the  iono¬ 
sphere  are  actively  contributing  to  the  signal  reception 
at  any  given  time  is  unknown.  During  sunrise  or  sunset 
periods  Yeh  and  Villard  ^  have  noticed  a  sinuous  type 
of  fading  which  they  attribute  to  the  simultaneous  recep¬ 
tion  of  the  signal  via  two  distinct  modes  of  propagation. 
It  has  been  speculated  that  one  of  the  modes  of  propaga¬ 
tion  is  caused  by  an  ionospheric  tilt  which  moves  west¬ 
ward  with  the  sun.  The  motion  of  the  tilt  imparts  a 
relatively  large  Doppler  frequency  shift  to  one  of  the 
received  waves.  The  other  wave  is  supposed  to  be  due 
to  ordinary  reflection  by  the  ionosphere.  The  two  waves 
then  combine  to  produce  a  single  v/ave  which  has  a  sinuous 

type  of  modulation.  This  effect  also  has  been  noticed  by 

(2) 

Nicholson. 

Pertinent  to  the  subject  of  fading  of  reflected  waves 
is  the  information  gained  by  Wild  and  Roberts  when 

studying  the  variations  of  amplitude  of  extra-terrestrial 
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noise  that  had  been  received  after  passing  through  the 
ionosphere.  In  this  case  the  ionosphere  is  assumed  to 
act  as  a  diffracting  screen  which  is  continually  chang¬ 
ing  its  structure  as  well  as  moving  around.  These  ex¬ 
perimenters  recorded  the  amplitude  of  the  received  signal 
versus  time  and  frequency,  using  a  radio  star  in  the  con¬ 
stellation  of  Cygnus  as  the  source  of  radiation.  The  re¬ 
ceiver  was  swept  over  the  frequency  range  of  40  to  70  me. 
The  period  of  time  required  to  sweep  the  receiver  was  short 
with  respect  to  the  time  required  for  changes  in  amplitude 
of  the  received  signal  at  any  given  frequency  in  the  range. 

From  these  measurements  they  were  able  to  determine 
that  the  amplitude  fading  appeared  to  be  created,  at  least 
some  of  the  time,  by  a  single  diffracting  irregularity  of 
large  dimension  which  was  moving  past  the  receiver.  The 
difference  between  the  conditions  of  this  test  and  the 
conditions  of  an  experiment  measuring  the  fading  of  a  re¬ 
flected  signal  are  many.  In  the  first  place  the  area  of 
the  ionosphere  involved  in  the  experiment  is  probably 
smaller  for  the  diffraction  type  of  experiment  than  for 
the  reflection  type  of  experiment.  Secondly,  there  is 
the  possibility  of  several  modes  of  propagation  existing 
simultaneously  in  the  normal  HF  radio  link.  However,  some 
results  that  are  useful  in  describing  the  fading  of 
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reflected  waves  have  been  obtained  in  the  diffraction  type 
of  experiment* 

The  motion  of  the  ionosphere  has  been  investigated  by 

( 4) 

a  number  of  researchers.  Briggs  and  Spencer  have 

mentioned  speeds  on  the  order  of  80  meters  per  second, 
using  the  method  of  spaced  receivers  which  was  developed 
by  S.  N.  Mitra.  There  is  some  doubt  as  to  the  accuracy 

of  these  figures  because  of  the  assumptions  inherent  in  the 
method  concerning  the  shape  of  the  irregularities.  Mot-ions 
measured  by  reflection  from  meteor  trails  establish  speeds 
on  the  order  of  70  meters  per  second  according  to  Manning 
and  Eshleman,  Measurements  o f  sodium  vapor  clouds  indi- 

cate  horizontal  wind  velocities  of  140  meters  per  second.  '' 
Velocities  reported  by  Wild  and  Roberts  derived  from  radio¬ 
star  scintillation  measurements  are  20  meters  per  second 
after  removal  of  apparent  motion  due  to  rotation  of  the 
earth.  ^  ^ 

The  time  of  growth  and  decay  of  small  scale  (30  km) 
irregularities  as  measured  by  a  backscatter  experiment  haB 

/  O  \ 

been  reported  by  Barry  to  be  ten  minutes  on  the  average. 

(9) 

Spencer  indicates  that  irregularities  measured  by 

virtue  of  their  diffraction  of  extra-terrestrial  waves  at 
38  me  were  elongated  by  a  factor  of  at  least  5:1  possibly 
along  the  direction  of  the  geomagnetic  field  lines. 
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The  motion  of  the  ionosphere  appears  to  bet  in  the 
large  sense,  a  uniform  motion  of  the  whole  reflecting 
portion  of  the  ionosphere.  However,  Court  reports 

that  there  is  some  randomness  in  the  varying  radio  waves 
that  cannot  be  explained  on  the  basis  of  a  uniformly 
translating  ionosphere. 

These  results  seem  to  suggest  that  only  a  small  num¬ 
ber  of  scatterers  is  usually  effective  in  causing  the  ampli¬ 
tude  of  reflected  waves  to  vary.  The  scale  of  turbulence 
seems  to  suggest  that  at  least  part  of  the  fading  is  caused 
by  an  interference  pattern  moving  past  the  receiver  that  is 
due  to  the  scattering  from  a  large  region  of  the  ionosphere 
moving  through  the  illuminated  area  of  the  ionosphere. 

4.  Mathematical  Models 

In  considering  mathematical  models  which  might  possibly 
describe  in  a  statistical  sense  the  varying  amplitude  of  the 
received  wave  we  might  consider  the  wave  received  at  the  re¬ 
ceiver  as  a  discrete  sum  of  received  signals,  each  signal 
being  thought  of  as  having  been  scattered  from  a  particular 
spot  whose  dimensions  are  all  smaller  than  a  wavelength.  In 
this  case  the  randomness  of  the  variation  of  the  amplitude 
of  the  received  signal  would  arise  from  variations  in  (a)  the 
number  of  such  point  scatterers  inside  the  intersection  of 
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the  transmitting  and  receiving  antenna  cones,  (b)  the 
strength  of  the  scattered  signal  from  each  of  the 
scatterers,  (c)  the  component  of  velocity  of  each  spot 
in  the  direction  of  transmission,  and  (d)  the  instan¬ 
taneous  phase  angle  of  the  return  from  the  individual 
spot. 


Another  possible  way  of  viewing  this  problem  is  to 
consider  a  Markov  chain  in  which  states  of  the  chain  are 
constructed  from  the  possible  diffraction  patterns  associated 
with  scattering  from  irregularities  which  have  dimensions  on 
the  order  of  many  wavelengths.  The  states  would  be  defined 
by  the  existence  of  any  one  of  these  patterns  or  of  combina¬ 
tions  of  these  patterns.  A  factor  which  must  enter  into  the 
delineation  of  the  states  is  the  speed  with  which  the  dif¬ 
fraction  patterns  are  moving  past  the  receiver.  This  type 
of  Markov  chain  is  too  complex  to  handle  in  a  simple  fashion. 

If  we  may  assume  that  the  process  is  stationary  we  can 
learn  something  from  the  power  spectral  density  of  the  ampli¬ 
tude  records.  If  there  is  no  power  density  above  say  10 
cycles  per  second  one  might  intuitively  expect  the  ionosphere 
to  be  relatively  constant  in  a  period  of  1/10  of  a  second. 

If  the  power  density  were  spread  uniformly  over  a  wide  range 
of  frequencies  one  might  expect  the  signal  was  composed  of 
many  random  components.  If  there  were  two  signals  adding 
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to  form  the  received  signal  such  that  one  of  these  compo¬ 
nent  signals  had  suffered  a  Doppler  frequency  shift  one 
might  expect  a  large  amount  of  fading  power  in  a  Gmail 
region  of  frequency  about  the  value  of  the  Doppler  fre¬ 
quency  shift. 

As  shown  later,  power  spectral  density  estimates 
made  from  fade  records  fluctuate  wildly  as  functions  of 
frequency.  The  ratio  of  a  maximum  value  to  the  adjacent 
minimum  value  may  be  as  high  as  three  for  a  typical  record. 
One  wonders  whether  these  peaks  in  the  spectral  estimates 
are  indicative  of  randomly  occurring  periodicities  in  the 
fade  records  or  whether  these  fluctuations  normally  occur 
as  deviations  of  the  estimates  from  the  true  value  of  the 
spectrum.  This  question  is  inextricably  related  to  the 
question  of  whether  adjacent  (in  time)  data  points  are 
statistically  independent  samples. 

The  sample  theorem  tells  us  that  when  the  true  power 
spectral  density  function  is  essentially  zero  above  some 
frequency,  f^,  the  time  period  between  samples  need  not  be 
less  than  2/f  .  However,  samples  taken  this  close  together 
are  obviously  not  independent  samples  but  are  very  well 
correlated.  On  the  other  hand,  if  we  take  samples  separated 
by  a  time  period  At  which  is  much  larger  than  2/f  so  that 
the  samples  are  independent  then  we  can  derive  no  information 
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concerning  the  shape  of  the  spectrum  between  the  fre¬ 
quency  f  =  l/(2At)  and  the  frequency 

When  the  sample  points  are  well  correlated  the  task 
of  finding  the  expected  value  and  the  standard  deviation 
of  the  spectral  estimates  (i.e.,  the  task  of  deciding  how 
greatly  the  spectral  estimates  might  vary  from  tne  true 
spectrum)  becomes  difficult. 

Despite  these  difficulties  the  time  increment  between 
sample  points  was  chosen  small,  enough  so  that  no  extreme 
amplitude  change  occurred  between  samples. 

5.  Digital  Computer  Program 

The  problem  of  estimating  the  power  spectral  density 
of  the  process  from  the  amplitude  records  was  handled  in 
two  different  ways.  One  method  was  to  simulate  a  low  fre¬ 
quency  spectrograph  of  adjustable  Q  with  an  analog  computer. 
The  other  method  involved  using  the  digital  computer  and 
discrete  data  points  to  estimate  the  auto- correlation  func¬ 
tion  and  then,  from  this  function,  to  estimate  the  power 
spectral  density. 

The  second  method  mentioned  was  used  to  obtain  only  a 
few  estimates  of  the  power  spectral  density.  The  program 
for  obtaining  these  estimates  is  described  in  a  Boeing 
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report.  The  estimates  of  the  auto-correlation  func¬ 

tion,  R(m),  were  computed  by  using  the  formula 


N-m 

R(a>  =  ^  H(nAt)  H  [(n+m)  At],  for  a  =,  0,  1,  2,..*,[N/4],  5.1 

n=l 

I  M 

where  [H(nAt)]  is  equal  to  the  difference  between  the 
data  (as  a  function  of  time)  and  the  average  value  of  the 
data,  and  where  N  is  the  total  number  of  data  points.  The 
estimates  of  the  power  spectral  density  S^(2nkAg)  were  then 
computed  by  the  use  of  the  formula 


[  N/4  ] 

Se(2nkAg)  =  4At  ^  R(m)  cos  (2n  kAg  mAt)J  , 

m=:l 

where  k  is  an  integer  and  Ag  is  the  increment  of  frequency. 
The  subscript  e  denotes  the  estimated  value  of  S.  These 
formulas  are  approximations  to  the  formulas 


r  t/2 

R(t)  =  1/T  f  (t)  f,  (t  +  t)  dt 

J  ~t/p  1  1 


and 


S(«) 


T/4 

0 


5.2 


5.3 


R(t)  cos  ( id t  )  dt 


5.4 
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No  weighting  function  wars  used  in  the  calculation  of 

the  auto-correlation  estimate  (i.e.,  in  the  terminology  of 

(12) 

Blackman  and  Tukey  ,  the  Tlag  window*  was  a  square  win¬ 

dow);  hence,  power  density  at  a  discrete  frequency  will 
appear  in  the  estimate  as  a  (sin  u)/u  function  centered 
about  the  frequency  of  the  power  density. 

The  digital  computer  program  did  not  turn  out  to  be 
as  useful  as  had  been  anticipated  since  the  amount  of  com¬ 
puter  time  required  to  handle  reasonably  long  data  samples 
was  much  too  great. 

6.  Statistical  Test  of  the  Spectral  Estimates 

Under  certain  conditions  one  may  derive  some  quantita¬ 
tive  information  concerning  the  true  power  spectral  density 
from  the  spectral  estimates  using  a  statistical  teBt  of  the 
magnitude  of  the  estimate.  That  is,  if  the  random  process, 
of  which  the  data  constitute  a  sample  function,  is  a  sta¬ 
tionary,  ergodic,  random  process  and  if  the  original  pro¬ 
cess  is  a  Gaussian  process,  then  a  probability  distribu¬ 
tion  defined  over  the  amplitude  of  the  power  spectral  density 
estimate  may  be  calculated.  From  this  probability  distribu¬ 
tion  one  may  find  the  probability  of  obtaining  any  given 
value  when  computing  the  power  spectral  density  estimate.  If 
the  spectral  estimate  obtained  has  a  small  probability  one 
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makes  the  inference  that  either  the  conditions  of  station- 
arity  and  ergodicity  or  the  requirement  that  the  process  be 
a  Gaussian  process  did  not  hold. 

The  test  described  above  is  modified  from  a  standard 
test  which  is  usually  applied  to  processes  in  order  to  dis¬ 
cover  periodicities.  As  a  matter  of  fact  the  standard  test 
does  not  necessarily  yield  information  about  the  frequencies 
or  existence  of  periodicities. 


The  statistical  reasoning  underlying  this  test  is 

described  by  the  following  argument.  Assume  that  the  data 

2 

are  samples  of  a  random  variable  with  variance  o  whose  sam¬ 
ple  functions  are  functions  of  time  X(t).  We  are  considering 
that  we  might  have  obtained  any  one  of  an  infinite  number  of 
possible  sample  functions  and  that  the  possible  sample  values 
X(nAt),  denoted  X^,  and  X(  nAt  +  mAt),  denoted  %n+mi  have 
the  joint  probability  density  function 


f(X  .  X  ) 
n*  n+ra 


2pXn  W  *2n+J\  ) 
2o2(l  -  p2)  )  ) 


6.1 


where  the  correlation  between  X  and  X  is  p. 

n  n+ra 

In  this  formulation  it  may  be  proved  that  the  quantity, 

p,  is  an  exponential  function  of  the  absolute  value  of  the 

(15) 

time  difference,  mAt.  It  is  also  possible  to  prove 

that  the  joint  probability  density  of  any  k  sample  values  is 
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multivariate  normal* 


In  writing  the  density  function  as  above  we  have 
implicitly  assumed  that  the  random  process  has  zero 
expectation,  a  condition  which  can  be  approximated  by 
subtracting  the  average  value  of  the  signal. 


If  the  total  length  of  the  record  is  denoted  by  T 
and  if  At  =  T/N  we  may  calculate  the  spectral  density 
function  of  a  given  sample  record  as  a  function  of 
radian  frequency  by  the  formula 


N 


2 


s  U,  )  =  *nrr 

e  k  NAt 


Z  Xn  exp  {“  Jwk  nAt}  At 
n=l 


6.2 


As  may  be  shown  this  is  equal  to 

N 

Se(w^)  =  4At  ^  R(m)  cos  (w^mAt)  ]  6.3 

m=l 

which  is  related  to  the  function  computed  on  the  digital 
computer. 

V/e  see  that  with  different  sample  records  we  will 
obtain  different  functions  se(w^)  B0  that  for  any  parti¬ 
cular  value  of  radian  frequency  we  can  think  of  as 
being  a  random  variable  with  some  probability  distribu¬ 
tion  defined  over  the  possible  values  of  8  . 


This  random 
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variable  is  related  to  the  set  of  random  variables  lx  V 

i nf 

where  n  ranges  from  1  to  N  so  that  the  probability 
distribution  for  is  related  to  the  probability  dis¬ 
tribution  for  X^.  That  is,  in  order  to  find  the  pro¬ 
bability  distribution  for  we  must  know  the  probability 
distribution  for  X^.  This  is  the  reason  for  assuming 
that  the  process  X(t)  is  a  Gaussian  process. 


As  may  be  shown 


(14) 


the  probability  distribution 


for  S  is  given  by 


f(S  ) 
e 


exp 


2  a 


{-=*!• 


p 

where  o'  is  the  variance  of  the  original  process  X(t). 
With  this  probability  distribution  we  may  compute  how 
likely  any  particular  value  of  Sg  is.  In  fact  we  may 
compute  the  probability  that,  for  a  particular  value  of 
radian  frequency,  say  is  less  than  some  number  C. 

This  probability  is 


6.4 


/ 

l  =  JLi 

fC 

(  NS  ) 

Jo  S 

exp 

'  0 

l  e 

’  2o2J 

l  2o2  } 

Nr 

dS  =  1  -  exp  (-  — . 

e  ~  2 

20 


6.5 


V/ hen  analyzing  the  spectral  estimates  for  a  given  record 
v/e  are  most  interested  in  peak  values  of  these  estimates 
(i.e.,  the  maximum  value  of  the  function  over  all 

values  of  u>).  We  would  like  to  know  what  significance 
this  peak  has.  In  other  words,  we  would  like  to  calculate 
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a  value  C  for  which  the  following  statement  is  true; 

p|max*  over  of  <  0.01  . 

This  may  be  accomplished  by  computing  the  value 
of  C  for  which 


6.6 


fSe(wl>  =  0  fSe(V  =  C  r  i 

iW .  Wj 

J  se(Wl)  =  o  J  se(»M)  =  0 


dS  (w.,)  ...  dS  (»J  =  0.99  6.7 

el  e  M 

where  . ..,  is  joint  probability 

density  function  for  the  M  values  of  that  have  been 
estimated. 

Tills  joint  probability  density  function  is  diffi¬ 
cult  to  derive  from  the  original  probability  density 

functions  for  X  and  X  because  of  the  fact  that 
v  n+ra 

and  are  not  necessarily  independent  ran¬ 

dom  variables  because  the  random  variables  X  and  X 

n  n+m 

are  not  necessarily  independent.  However,  for  a  set  of 
carefully  selected  values  of  k  it  is  theoretically  pos¬ 
sible  to  show  that  the  random  variables  S  (w,  )  for  k  in 

6  K 

this  set  are  independent  random  variables.  In  that  case 
the  joint  probability  density  function  required  in  the 
integral  above  factors  into  the  product  oC  the  separate 
probability  density  functions,  and  the  integral  factors 
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into  the  product  of  M  identical  integrals. 

Finally  the  desired  probability  is 

P  |  max  over  of  >  c|  =  1  -  -  exp  (-  -~~f)  J  <  0.01.  6.8 

This  final  relation  may  be  solved  for  the  desired  value  of 
C  if  the  variance  of  the  original  process  is  known. 

The  test  just  described  is  a  slight  modification 
of  the  standard  tests  (see  M.  G.  Kendall  differing 

because  there  is  no  assumption  made  as  to  the  independence 
of  the  original  sample  points.  As  indicated  previously, 
if  the  sample  points  are  widely  spaced  so  that  the  sample 
points  are  statistically  independent  then  we  cannot  esti¬ 
mate  the  spectral  density  function  for  frequencies  suf¬ 
ficiently  high  as  to  be  of  great  interest. 

The  test  described  above  also  differs  from  one  des¬ 
cribed  by  Grenander  and  Rosenblatt  because  the 

selected  frequencies  which  may  be  used  are  different. 

Finally  the  test  described  above  has  the  drawbacks 
that  the  variance  of  the  original  process  must  be  esti¬ 
mated  (with  unknown  probability)  and  the  expected  value 
of  the  original  process  must  be  estimated  (again  with 
unknown  probability)  before  the  value  of  the  confidence 
level  may  be  computed. 
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On  the  credit  side  the  test  described  above  is 
independent  of  the  correlation  of  the  original  sample 
points  and  hence,  if  the  confidence  level  is  exceeded 
we  may  infer  that  either  the  original  process  was  not 
stationary  and  ergo die  or  that  it  was  not  Gaussian  or 
both  without  having  to  include  also  the  possibility  that 
the  original  sample  points  were  not  independent. 


7 .  Statistical  Properties  of  the  Spectral  Estimates 

We  need  to  investigate  the  statistical  properties 
of  the  random  variables  S ^ ) .  To  do  this  we  start  by 
examining  the  properties  of  two  functions 

N 


and 


A(cok)  r=  ^  X(nAt)  cos  (w^nAt)  At 


n=l 


nr 

B(o)^)  ^  X(nAt)  sin  (w^nAt)  At 


n=l 


which  are  related  to  S  (or  )  as  below: 

e  k 


•W  '  A?(V  *  Bp("k)* 


We  should  determine  those  values  of  uj.  for  which 

k 

Se(w^)  are  a  set  of  independent  random  variables. 
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We  first  note  that  ACu^)  and  B(oj^)  are  composed  of 

a  linear  combination  of  the  random  variables  X  .  Since 

n 

we  have  assumed  that  the  random  variables  X  are  Gaussian 

n 

we  also  have  ACu}^)  and  B(w^)  as  Gaussian  random  variables. 

Further,  because  they  are  Gaussian,  we  can  establish  the 

independence  of  ACu^)  an(l  A(oj^)  or  simply  by  showing 

that  the  correlation  between  these  random  variables  is 

zero.  If  the  two  random  variables  A(oj^)  and  BCtu^)  are 

independent  of  the  two  random  variables  A(cu^)  and  B(to^) 

then  the  random  variables  S  (w.  )  and  S  (w.)  must  be  in- 

e  k  e  *• 

dependent . 

An  incidental  result  we  should  obtain  in  establishing 
the  allowed  values  of  is  an  expression  for  the  expecta¬ 
tion  of  S  (w,  )  as  a  function  of  u,  It  is  known  that  the 
e  k 

expected  value  of  the  estimate  approaches  the  true  power 
spectral  density  as  N  approaches  infinity.  However, 

the  expected  value  for  finite  N  would  give  us  some  idea  of 
what  variations  we  might  expect  in  the  spectral  estimates 
over  the  entire  set  of  frequencies. 

We  see  that  E^S^Co^) J  =  E^A£*(w^)J  +  and 

also  we  see  that  the  correlation  between  ACoj^)  and  A(to^) 
is  given  by 

E[A(uk)  AU£)]  -  Ej  A(wk)]E[A(wx)] 

i/E|V(tok)  e[a2(o)£) 


7.4 
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so  for  both  objectives  we  must  compute  the  covariance, 
0^,  of  the  random  variables  A(u>k)t  A(w^),  B(wfc)  and 
B(w^) . 


Th 


e  expectation  of  the  product  of  A(u)^)  and  A(w^) 


is 


E^A(u)k)  A(o)^)J  =  E  ["  ^  X(nAt)  cos  (a^nAt)] 

n=l 

N  i 

["^X(jAt)  cos  (co^jAt)jj 


j=l 


N 


^  E^X(nAt)  X  (jAt)J  cos  (to^nAt)  cos  (w^jAt), 
ni  J=1 


7.5 


If  the  process  is  ergo  die  and  if  E[X(nAt)]  is  zero  we 
have  E[x(nAt)  X(jAt)]  aB  the  covariance,  and  also 

the  equality 

°kx=  E[A(V  a<«p] 

2  N 

=  V  p(|n-j|At)  cob  (w^iiAt)  cos  (u^jAt)  7.6 

n,  j=l 

2 

where  0  is  the  variance  of  the  X(t)  process  and  p  is 
the  correlation  between  X(nAt)  and  X(jAt).  As  stated 
before,  if  the  process  is  a  Markov  Gaussian  process  it 
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can  be  proved  that  p  has  the  form 

/  I n- j | At  ) 
p  -  exp  <  -  J ^ j 

where  t  is  a  constant  dictated  by  how  wide  a  frequency 
band  the  original  process  occupies. 


By  rearranging  the  order  of  summation  we  can  obtain 
^  N-l  N 


_  4o_ 

Jk i  “  NAt 


^  ^  p(uAt)  cos  (w^XAt)  cos  j^ca^CX-u)  AtJ  7.^ 


u=l  X=u+1 

and  then  by  the  use  of  trigonometric  identities  we  have 
N-l  N 


2  1  / 

ak,e  =ftt  X  I  p(uAt)  |COB[x(“k  +  "s)  ^  -  “V^tJ 

u=l  X=u+1  " 

+  cos  -  w£)  XAt  +  w^uAtjj 


7.9 


which  may  be  evaluated  to  give 


N-l 


“  NAt 


u=l 


r ( N  +  1)  rr  +  upj  sin  j^(N  -  u)  aj 


^  p(uAt) 

1 

cos  j^(N  +  l)p  +  uaj  sin  |*(N  -  u)pj 


where 


sin  P 


(“k+  “P  ..  _  <<v  <■»*> 

a  =  - 5 -  At,  p  =  - - -  At. 


7.10 


7.11 


If  p(t)  is  known  then  the  last  sura  can  be  evaluated 
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hut  it  is  advantageous  to  examine  the  terms  as  they 
stand.  Upon  examination  of  the  auto-correlation  esti¬ 
mates,  Figures  3  and  4,  it  appears  probable  that  the 
true  correlation  drops  to  zero  in  a  time  less  than  one 
second;  hence,  x  is  lees  than  one.  This  means  approxi¬ 
mately  the  first  five  terms  in  the  sum  are  the  signifi¬ 
cant  terms. 


We  cannot  expect  to  estimate  the  spectrum  for  fre¬ 
quencies  higher  than  the  frequency  for  which  we  have  two 
sample  points  per  cycle.  This  means  that  the  maximum 
value  of  a  was  chosen  to  be 


Max. (w^At ) 


2n  At  ^  2 re  At 

T,  -  2  At  = 

k 


so  that  a  lies  between  re  and  zero.  In  fact  we  have 
chosen  our  sampling  time  so  that  a  lies  between  u/2  and 
zero . 


7.12 


As  a  approaches  zero  we  see  that  the  sum  approaches 


small  u 


^  p(uAt)  | 


cos( Nq )  sirt( No) 
sin  a 


cos(Np ) sin(Nft) 
sin  p 


u=a 


7.13 


or 


small  u 

Y  p(«At)  1/?  {giSlgWgQ 
L  {  sin  a 

u=l 


7.14 
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Consequently  we  would  expect  the  cross  correlation  to 
fluctuate  as  a  function  of  a  and  3  between  zero  and 

some  large  number  whose  magnitude  is  less  than 

p small  u 
2o  V 

— )  p(uAt)  and  whose  magnitude  varies  as 
u  =1 

|  l/sin(a)|  +  |l/sin  (p)|  as  Na  and  N8  pass  through 
the  points  M  n  and  (2M  +  l)n/2  (sin  Na  =  0  and 
sin  Na  =  1). 


The  expected  value  of  depends  upon  and 

we  see  that  is  given  by  the  Eq.  7.10  where  p  is  zero 

and  a  is  w^At.  In  that  case  we  have 


kk“ 


2o 

NAt 


small  u 
^  p(uAt) 
u=l 


|^(N+l)aJ  sin  |^(N~u)aJ 
sin  a  . . . 


+  cos  (ua) 


7.13 


and  it  is  apparent  that  this  fluctuates  wildly  as  a  func¬ 
tion  of  a  =  w^At  ^or  sma-^  a* 


To  find  the  allowed  values  of  w.  the  summation  should 

k 

be  performed  and  the  resulting  transcendental  expression 
should  be  set  to  zero.  Then  the  resulting  equation  should 
be  solved  for  the  allowed  values,  . 

One  of  the  conclusions  to  be  drawn  from  this  analysis 
is  that  the  value  of  the  minimum  frequency  at  which  the 
spectral  density  is  estimated  should  be  greater  than  the 
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frequency  for  which  about  five  or  ten  degrees. 

If  the  process  contains  periodicities  (as  might 
be  the  case  when  two  modes  of  propagation  exist  sim¬ 
ultaneously)  we  would  expect  peaks  occurring  in  the 
spectral  estimates  at  the  frequencies  of  the  periodi¬ 
cities.  The  width  of  these  peaks  would  be  related  to 
the  duration  of  the  periodic  component  or  to  the  dura¬ 
tion  of  the  record,  whichever  was  shorter.  The  purpose 
behind  not  smoothing  the  spectral  estimates  by  the  use 
of  a  weighting  function  in  the  data  was  to  produce  the 
narrowest  possible  peaks  from  these  periodicities. 

Suppose  that  the  record  lasted  for  T  seconds  and 
that  it  contained  a  sinusoidal  component  that  lasted 
for  L  seconds  where  T  is  the  order  of  five  or  ten  times 
L  and  where  there  are  not  less  than  about  five  cycles 
of  sinusoidal  fade  in  evidence.  We  then  can  charac¬ 
terize  the  received  signal  amplitude  as  being  the  sum 
of  a  sinusoidal  function  and  a  random  function,  X(t), 
as  follows: 

f(t)  =  A  sin  (at)  U(t)  U(L  -  t)  +  X(t) 
where  U(t)  is  a  unit  step  function.  In  that  case  i * 


7.16 


can  be  shown  that  the  expected  value  of  the  spectral 


estimate  is  approximately 


fT/4 

+  p(r)  cos  (no)  t)  dx  7. 

Jo  U 

where  the  f i ret  term  Is  due  to  the  sinusoid*  This  ex¬ 
pression  displays  the  peaks  mentioned  above. 

The  identification  of  these  peaks  remains  a  problem 
whose  solution  depends  upon  a  more  complete  knowledge  of 
the  possible  variations  of  the  spectral  estimates.  Know¬ 
ledge  of  the  values  of  w  at  which  the  cross- correlation 
between  spectral  estimates  is  zero  would  help  considerably; 
however,  the  worth  of  pursuing  this  line  of  work  is  ques¬ 
tionable. 

The  digital  computations  are  more  accurate  than  the 
analog  computations,  but  this  appears  to  be  the  only  ad¬ 
vantage  the  digital  computer  program  has  over  the  analog 
computer  program.  Unfortunately  the  length  of  time  re¬ 
quired  to  compute  spectral  density  estimates  with  the 
digital  computer  increases  approximately  with  the  square 
of  the  number  of  data  points*  The  five  minute  sample  of 
data  used  for  the  estimate  shown  in  Figure  8  (page  4?) 
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contained  600  points  and  required  4.6  minutes  of  computer 
time.  The  full  22  minutes  of  data  in  this  particular  sam¬ 
ple  would  require  about  an  hour  and  a  half  of  computer 
time.  At  this  rate  the  IBM  704  computer  would  require 
about  9*3  hours  to  analyze  a  55  minute  sample  such  as  has 
been  analyzed  with  the  analog  computer  in  about  one  hour. 

Because  the  analog  computer  can  handle  longer  sam¬ 
ples  of  data  the  power  spectral  density  estimates  are 
considerable  more  stable;  that  is,  they  donft  fluctuate 
as  Eq.  (7.15)  predicts  that  the  digital  computer  esti¬ 
mates  will  fluctuate.  This  fluctuation  is  shown  in 
Figure  8  where  the  digital  computer  estimate  of  the  power 
spectral  density  for  data  in  the  time  interval  of  1431 :l8 
to  1436:20  PST  on  August  4,  i960  is  graphed  along  with 
the  analog  computer  estimate  of  the  power  spectral  density 
for  the  data  in  the  time  period  from  1421  to  1443  PST  on 
August  4,  i960. 

The  relative  calibration  between  the  amplitudes  of 
the  two  estimates  is  unknown  so  only  the  general  shape 
of  the  curves  can  be  compared.  When  comparing  the  two 
curves  one  must  keep  in  mind  that  the  digital  computer 
estimate  is  for  five  minutes  of  the  data  whereas  the 
analog  computer  estimate  is  for  the  whole  22  minutes. 
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8 .  Analog  Computer  Program 

The  purpose  of  the  analog  computer  program  is  to 
obtain  spectral  estimates  from  the  sample  functions  and 
to  identify  periodic  fades  and  to  obtain  the  following 
information  concerning  these  periodic  fades:  (1)  the 
duration  time,  (2)  the  starting  and  ending  time,  (3) 
their  variation  in  ‘frequency*  and  (4)  their  amplitude. 

To  obtain  this  information  from  a  known  sample 
function,  e^Ct),  the  analog  computer  is  programmed  to 
solve  the  second-order  equation 

e!,+  2  o  e‘  +  cr^)  e  =  <u^e  8.1 

o  O  0  3 

where  w  and  a  are  constants  that  are  adjusted  by 
changing  a  dc  voltage  in  the  computer.  The  voltage, 
eQ(t),  is  then  similar  to  the  output  of  a  tuned  circuit 
that  has  been  excited  by  the  voltage  e^(t).  The  average 
of  the  square  of  the  voltage,  eQ(t),  is  computed  as  a 
measure  of  the  power  density  in  e  (t)  at  the  frequency 

•3 

w.  A  diagram  of  the  computer  circuit  1r  shown  in  Figure 

3. 

The  data  are  treated  in  the  following  way.  The 
signal  is  recorded  on  magnetic  tape  at  0.073  inches 
per  second  and  then  played  back  at  7.3  inches  per  second 


Analog  computer  circuit’ 
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and  re-recorded  on  a  25  foot  loop  of  magnetic  tape.  Thus, 
the  time  scale  is  compressed  by  a  factor  of  100.  Along 
with  the  data  a  synchronizing  signal  which  occurs  at  the 
same  time  as  the  start  of  the  data  is  recorded  on  an  ad¬ 
jacent  channel  of  the  magnetic  tape.  This  tape  is  then 
reproduced  with  a  loop-tape  magazine  so  that  the  data  are 
reproduced  over  and  over  again  in  a  repetitive  fashion* 

On  the  first  reproduction  of  the  data  the  voltage  out 
of  the  magnetic  tape  recorder  is  squared  in  the  analog  com¬ 
puter  and  is  put  into  the  simulated  resonant  circuit  of  the 
computer*  The  simulated  resonant  circuit  is  tuned  to  a  low 
frequency  of,  say,  one  cps.  The  synchronizing  signal  starts 
a  pre-set  timing  circuit  in  the  computer.  Samples  of  the 
computer  outputs  are  shown  in  Figure  6. 

At  the  end  of  the  first  reproduction  of  the  data 
the  computer  de-energizes  its  circuits,  reads  out  the  aver¬ 
age  of  the  square  of  the  response  of  the  simulated  tuned 
circuit  and  resets  the  tuning  of  the  tuned  circuit  to  a 
higher  frequency,  say,  1.03  cps.  At  the  start  of  the  second 
reproduction  of  the  data  the  synchronizing  signal  re-act i- 
vates  the  computer  circuits  and  the  cycle  is  repeated.  In 
this  fashion  a  measure  of  the  power  density  is  obtained  for 
frequencies  spaced  throughout  the  range  from  0.01  cps  to 
1  cps  (as  reckoned  on  the  time  scale  of  the  original  data)* 
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Fig.  6.  Analog  computer  output  voltages0 
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In  order  to  briefly  outline  the  properties  of  the 
spectral  estimates  obtained  from  the  analog  computer  we 
may  analyze  the  differential  equation  which  the  computer 
solves.  This  equation  is 

e,!+  2a  e'  +  (a>^  +  C^)  e  =  e  8.2 

O  O  O  o  o  Ci 


where  e  is  the  square  of  the  data  waveform  and  e  is  a 
e  o 

voltage  yielding  a  measure  of  the  power  density.  If  eg 
were  a  sinusoid  of  frequency  co^  the  voltage  would  ap¬ 
proach  a  sinusoid.  We  may  find  this  sinusoid,  by  finding 
the  steady  state  solution  of  the  differential  equation  by 
solving  the  equation 

-on 2  E  +  2jo)^o  E  +  (w  o^)  E  =  w  ^  E  •  8.3 

L  O  U1  O  O  O  O  B 

o 

Solution  of  this  equation  gives  us 


E  =r 
O 


J,  2  2  2^2  n  2 

V  (u>  +  a  -  w-  )  +  C 2u>-.  a) 

oil 


8.4 


where  is  the  magnitude  of  the  sinusoidal  voltage  eQ. 
We  may  find  the  peak  value  of  the  response  when  the 
frequency  is  varied  by  taking  the  derivative  of  thiB  ex¬ 
pression  with  respect  to  u>^  and  then  by  setting  this  de¬ 
rivative  equal  to  zero.  The  result  of  this  is  that  the 

I  2  2 

response  Eq  is  maximum  at  vwQ  “  °  •  When  this  con¬ 

dition  is  met  the  formula  for  Eq  may  be  rewritten  as 
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E  = 
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v/2 


QE6 


v/2"QE6, 


8.5 


0) 

where  Q  =  •—  and  we  see  that  we  can  change  the  tuning  of 

2 

the  resonant  system  without  changing  the  ratio  of  E  /E 

o  s 

if  the  Q  is  kept  constant.  The  analog  computer  was  set 

so  that  Q  =  2ti  for  the  present  data.  With  this  or  higher 

values  of  Q  we  have  the  resonant  frequency  of  the  system 

f  p  p 

V -  o  equal  to  ui^  within  l.)l  The  change  in  tuning 
of  the  resonant  system  which  is  made  between  runs  of  the 
data  is  made  as  a  percentage  of  the  previous  resonant 
frequency.  The  frequency  increment  added  to  the  reso¬ 
nant  frequency  is  equal  to  one  half  of  the  bandwidth  of 
the  tuned  system.  Because  of  this  way  of  changing  fre¬ 
quency  the  resultant  plots  of  power  density  versus  fre¬ 
quency  are  made  on  logarithmic  graph  paper. 


If  eg  is  some  function  other  than  a  periodic  func¬ 
tion  we  may  find  an  integral  solution  of  the  differ¬ 
ential  equation  by  the  use  of  the  Laplace  transform. 
Solving  by  the  use  of  the  Laplace  transform  and  using 

the  information  that  e  (0)  and  e’(0)  are  both  equal  to 

o  o 

zero,  we  find  that 

a)  2  E  (s) 

E  (e)  =  — - 2 - 6.  —  . 

o  c  ^ 

s  +  2as  +  (a)  +  a  ) 

o 


8,6 


Inverting  this  expression  with  the  use  of  the  convolu¬ 
tion  integral,  we  obtain  an  integral  form  of  the  solu¬ 
tion 


e  ^ 
o 


2 


co 


o 


t 

eg(X)  exp  |^-  a  (t-X)J  Gin  jjtfQ(t-X)J  dX  . 
o 


This  is  an  exact  expression  for  the  computed  voltage 

during  the  time  the  computer  is  1  tuned.1  to  the  frequency 

co  . 
o 


9.  Expected  Value  of  the  Spectral  Estimates 

The  measure  of  power  density  contained  in  a  given 
record  at  a  frequency  u>  is  computed  in  the  computer  by 
the  following  formula 

T 

Se=  6o2(t)  dt- 

J  o 

We  use  this  number  as  an  estimate  of  the  power  density 
of  the  true  power  spectrum  of  the  random  process  which 
produced  the  data  record.  In  order  to  decide  whether 
this  is  a  good  estimate  of  the  value  of  the  true  power 
spectrum  we  may  compute  the  expected  value  of  and 
discover  how  the  power  density  of  any  given  record  is 
related  to  the  true  power  spectrum. 

In  finding  the  expectation  of  we  regard  the 

input  signal  e  (t)  as  a  random  variable  so  that  we  are 
s 
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averaging  over  all  inputs  that  we  might  possibly  get 
from  the  random  process  with  the  spectral  density 
F(io)*  The  desired  expected  value  is  obtained  from  the 
calculation 

T 

E  [^e]  =  wT  P  E[V(t)]dt-  9.2 

J  o 


Hence,  we  must  compute 
.t 


E  [e02(T)]  =  E  FI  eg(X)  exp  [-o(t-X)]  sin  [io(t-X)]  dxj 
f  eg(r)  exp  [-a(t-x)]  sin  [(o(t-x)]  dxj 

J  o 

Combining  this  into  a  double  integral  and  then  using 
trigonometric  identities  and  the  following  definition 
of  the  correlation  function 

E  fe  (X)e  (t)1  -E2[ e  (t)] 

P(x-x )  - | — =! - Ll - =L 


9.3 


9.4 


where  is  equal  to  E  j^eB2(t)j  ~  E2  an(*  where 

j” c? s ( t. )  i8  assumed  to  be  zero,  we  can  rewrite  the  ex¬ 
pected  value  as 


2  2  „t  t 


Po2(t)]=  2X  U.  p(X-x)  exp  [a(X+x~2t)] 


•  [cos  w(X-x)  -  cos  co(X+x~2t]  dxdX  . 


9.5 


The  variables  of  integration  in  the  above  double 
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So  that  we  may  relate  this  expression  to  the  power 
spectrum  of  the  random  process  we  substitute  for  the 
correlation  function  the  following  relation 


p(u) 


f 

*  o 


)  cos  (w^u) 


We  can  then  invert  the  order  of  integration  and 
evaluate  the  resulting  integrals  to  obtain 


E 


[eo2(t)] 


2  2 

co  a 

x 


exp  (-2at) 


J  o  a  +  (co., 


+co)  ‘ 


f  F(co  )  a  sin  (2cot)  +  (w1 +co) 

- 1 - dwl  - - 1 - 2 - 

+  (co-.-ai)  J  CO*  a  +  (to-,+0)) 

1  o  i  1 

pc*  l 

"  n  1 


cos  (2cot) 


dun 


F(co-,  )  a  sin  (2cot)  -  (io^-ui)  cos  (2a>t) 


c?2+  (co^-to)2 


(o2a  2  [  C°  a  sinKon +u>)t  -(un +co) cos  (an +co )t 

+  - *-  exp(-at)  F(u). )  - 1  dco 

2  (Jo  (co^+w)  j a  +  (to^co)  J 

+  |  f(W;l) 


a  sin 

-p 

3 

H 

3 

i _ I 

-  (w^-co)  COB 

(w  -u>)t  J 

a2+  (co^-co)2 

F(io 


/: 

i; 


F(u)1) 


a  sin  -  (w^+uj)  cos 


i  o  sin 

2  ,  v  2 

a  +  (u)^+w) 

^(u)^40>)tj  -  (tO-^-w)  COB 

|^(o)^+0))t 

o  sin  £ 

a  +  (to^-co) 

(o)^-uo)oj  +  (co^+to)  COB  £ 

(o)1+to)  t  j 

(co^+w)  j^a2+  (w1+co)2j 

dun 


dun 


dun 
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a  sin 

(w^-Uj)t 

+  ((0^-U))  cos 

(a)1-0))t  J 

ss}(s2i  (ay-o,)2’ 

/no 

o 


F(oj^)  0  sin  j^(w1+o))tj  +  (o)^+w)  cos  ^(w^+toHj 


o^+  (u^  +co )  ^ 


dm. 


r 00  F(w^)  0  sin  (u)^-a))tj  +  (ca^-co)  cos 

( <n .  -  w)  t 

-  1  .  .1  j.. 

f  2  (  x2 

J  O  0)^  0  +  (u)j-(l)) 

2  2^ 

a)  o  rm 
x 


/r  w/(i)1  to/u)  n 

F(u>i)  h  7  72  “  “2  7  72 J  dwi 

o  a  +  (to  -u>J  0  +  (u)  +w) 


In  this  expression  the  first  three  terms  have  the  factor 
exp  ( - 20 1 ) ,  the  next  eight  terms  have  the  factor 
exp  ( —  ot )  but  the  last  term  contains  no  such  factor.  The 
significance  of  this  fact  is  that  after  ten  cycles  of  the 
deeired  frequency  each  of  the  first  eleven  terms  is  multi¬ 
plied  by  a  factor  that  is  less  than  0.01.  When  the  re¬ 
cords  are  about  twenty  cycles  or  more  of  the  desired  fre¬ 
quency  in  length  we  may  approximate  the  expression  for 
the  expected  value  with  just  the  laBt  term  of  the  fore¬ 
going  equation.  Then  the  expected  value  of  the  spectral 
estimate  is 


9.10 


9.H 
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The  integrand  in  this  expression  is  a  product  of  the 
function  describing  the  power  spectrum,  F(to^  )t  and  a  weight¬ 
ing  function.  A  plot  of  the  weighting  function  is  shown  in 
Figure  7.  The  weighting  function  has  a  peak  at  the  fre¬ 
quency  at  which  the  spectrum  is  being  estimated  and  drops 
rapidly  to  zero  on  either  side  of  this  peak.  If  we  define 
the  side  frequencies  as  those  frequencies  at  which  this 
plot  is  down  by  1/2,  we  see  that  we  have  a  ’bandwidth’  of 
roughly  l/Q,  where  Q  is  the  figure  of  merit  of  the  simulated 
tuned  circuit. 

This  form  of  weighting  function  assures  that  the 
greatest  contribution  to  the  expected  value  of  the  esti¬ 
mate  comes  from  F(w^)  in  the  neighborhood  of  tho  frequency 
of  the  estimate. 

As  a  check  on  the  work  we  may  compute  the  form  of  the 
estimate  for  special  types  of  Ftco^).  If  the  true  power 
spectral  density  function  is  a  constant,  A,  over  the  whole 
spectrum  then  the  expected  value  of  the  estimate  is 


9.12 
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Then  we  may  change  variables  as  per  the  relations 
o)^/ to  —  u,  dto^  =  a)  du 

and  the  expected  value  of  the  estimate  becomes 

1  _  1 _ 

u-l)2+  -I5  (u+l)2+  -ir 

4<T  4q2 

Since  c 0  does  not  appear  in  the  integral  the  estimate  is 
a  constant  proportional  to  the  true  value  of  the  spectrum. 

10.  Results 

Two  samples  of  the  data  for  which  spectral  estimates 
have  been  made  are  shown  in  Figures  1  and  2.  The  first 
sample  in  Figure  1  is  a  recording  of  the  detected  voltage 
from  a  Collins  5U4  receiver  lasting  for  22  minutes  ob¬ 
tained  on  August  4,  i960.  The  receiver  was  adjusted  so 
that  the  detected  voltage  was  directly  proportional  to  the 
r.f.  signal  voltage.  The  signal  was  received  at  17*9  me 
over  a  6000  km  path  arriving  at  approximately  13  degrees 
south  of  east.  The  second  sample  in  Figure  2  is  a  record¬ 
ing  of  the  detected  voltage  from  the  IF  of  a  Collins  75S1 
receiver  lasting  for  56  minutes  obtained  on  May  9,  1962. 
Again  the  detected  voltage  is  proportional  to  the  input 
signal.  The  signal  (WWV)  was  received  at  20  me  over  a 
3710  km  path  arriving  at  an  angle  of  about  2  degrees  north 


E  (Se)  =  -  (Aox2/2)J  (1/u)[- 

J  o  L( 


9.13 


9.14 
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of  east.  The  direction  of  increasing  signal  strength  in 
Figure  1  is  downward  whereas  it  is  upward  in  Figure  2. 

In  both  samples  there  is  considerable  variation  in 
the  amplitude  of  the  received  signal.  The  signal  actually 
fell  below  the  noise  level  of  the  receiver  on  occasion. 
However,  the  rate  of  variation  is  not  excessive.  The 
fades  do  not  exceed  one  cycle  per  second  in  either  sample. 
In  a  few  cases  (not  analyzed)  fades  of  one  to  two  cycles 
per  second  have  been  observed. 

The  firet  data  sample  appears  to  contain  no  fades 
which  oould  be  described  as  being  caused  by  a  beat  between 
signals  arriving  via  two  modes  of  propagation.  The  second 
sample  has  such  fades,  one  occurring  between  0050  and  0035 
GMT.  There  appears  to  be  a  definite  beat  of  0.3  cycles 
per  second  lasting  for  tens  of  cycles.  Another  appears 
between  0039  and  004l  GMT  with  a  frequency  of  0,15  cps. 
Both  of  these  beats  appear  to  be  superimposed  onto  a 
slower  type  of  fade  which  has  larger  peak  values. 

The  data  in  Figure  2  contain  many  peaks  of  signal 
laBting  20  and  30  seconds.  Between  2355  and  235^  GMT 
there  appears  to  be  a  periodic  increase  of  signal  with 
a  frequency  of  approximately  0.025  cps.  On  the  other 
hand,  this  type  of  fade  is  lacking  in  the  data  shown  in 


Figure  1 
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Wild  and  Roberts^ ^  observed  that  some  of  their 
scintillation  data  appeared  to  be  caused  by  a  single  lens¬ 
like  irregularity  in  ionization  which  caused  a  single  large 
increase  in  the  signal  received.  It  is  possible  that  a 
similar  situation  accounts  for  the  longer  lasting  peaks  of 
signal  strength  in  the  present  data* 

The  spectral  density  estimates  computed  with  the 
IBM  704  digital  computer  are  shown  in  Figures  8  and  9. 

The  estimate  chov/n  in  Figure  8  is  for  a  five  minute  section 
of  the  data  shown  in  Figure  1*  The  time  interval  of  the 
section  analyzed  lies  between  1431:18  and  1436:20  PST*  The 
estimate  of  the  spectral  density  function  is  the  solid-line 
curve*  The  estimate  shown  in  Figure  9  is  for  data  taken 
on  June  15,  I960,  from  l6l6:09  to  l6l8  PST*  These  data  are 
a  recording  of  the  signal  voltage  out  of  the  detector  of  a 
R390A/URR  receiver  mentioned  in  the  equipment  appendix.  The 
signal  was  transmitted  over  a  6000  km  path  at  11.9  nic. 

Both  of  these  estimates  exhibit  the  large  fluctuation 
at  low  frequencies  which  tendB  to  drop  off  at  higher  fre¬ 
quencies,  say  0*2  cps.  If  these  fluctuations  are  caused 
solely  by  the  rectangular  truncation  of  the  data  they 
should  have  a  period  of  0*0131  cps  for  the  estimate  in 
Figure  8  and  a  period  of  0.038  cps  for  the  estimate  in 
Figure  9*  This  period  is  directly  related  to  the  time 
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Fig,  8. 


Power  spectral  density  estimate  of  the  fading 
power  in  the  detected  17.9  mo  signal,  August  4, 
I960,  as  computed  by  digital  and  analog  methods. 
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length  of  the  auto-correlation  function  computed  for  the 
records.  Since  this  is  computed  for  a  time  period  of 
one-fourth  of  the  length  of  the  record  (see  equation  (5-1)) 
we  can  compute  the  period  of  the  oscillations. 

There  are  two  peaks  in  the  estimate  of  Figure  8  which 
appear  to  he  wider  than  the  remaining  peaks,  one  at  about 
0.115  cps  and  one  at  about  0.26  cps.  These  peaks  have  a 
width  at  their  bases  of  approximately  0.04  cps  and  0.035 
cps  respectively.  If  these  peaks  were  caused  by  sinusoids 
(or  quasi-sinusoids)  we  may  measure  their  effective  width 
and  compute  a  value  for  the  quantity  L  in  equation  (7*17) 
where  L  Is  the  duration  time  of  the  sinusoid.  On  this 
basis  the  sinusoids  must  have  lasted  about  25  seconds  and 
30  seconds  respectively.  That  would  correspond  to  3  cycles 
and  8  cycles  of  the  respective  frequencies.  This  inter¬ 
pretation  of  these  peaks  is  not  satisfactory  as  one  can 
see  when  inspecting  the  raw  data  (Figure  1).  The  analog 
computer  estimates  are  considerable  more  useful  since  they 
do  not  fluctuate  as  wildly. 

The  spectral  estimate  shown  in  Figure  9  shows  power 
density  at  considerably  higher  frequencies  than  that  of 
Figure  8.  Unfortunately  the  data  for  this  sample  exiBt 
only  in  the  form  of  punched  cards  (not  on  magnetic  tape) 
and  hence,  cannot  be  readily  analyzed  on  the  analog  computer. 
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When  computing  the  spectral  density  estimates  with 
the  digital  computer  the  auto-correlation  function  is 
also  computed*  The  auto-correlation  functions  for  the 
data  of  August  4,  I960,  from  1431 :l8  to  1436:20  PST  and 
of  June  15j  I960,  from  l6l6:09  to  l6l8  PST  are  shown  in 
Figures  3  and  4.  These  plots  both  start  out  at  a  value 
equal  to  the  mean-Bquare  value  of  the  data  (zero  time 
lag).  The  curves  then  drop  rapidly  to  zero  and  there¬ 
after  oscillate  randomly  about  zero.  Data  points  sep¬ 
arated  in  time  by  a  period  greater  than  the  time 
required  for  the  auto-correlation  curve  to  drop  to  zero 
are  more  likely  to  be  statistically  independent.  The 
sustained  oscillations  are  apparently  indicative  of  a 
periodic  process. 

The  most  informative  and  useful  graphs  are  the 
power  spectral  density  estimates  computed  with  the 
analog  computer.  These  are  shov/n  in  Figures  10  through 
15.  All  except  one  of  these  estimates  start  dropping 
off  to  zero  in  the  neighborhood  of  0*2  cps  and  increase 
or  stay  constant  below  this  frequency,  showing  that  the 
power  density  of  the  fades  lies  below  0.2  cps.  Actually 
there  is  power  density  at  much  lower  frequencies,  say 
one  to  three  orders  of  magnitude  below  0.02  cps  since 
the  signal  fades  in  and  out  over  hour  periods  and  since 
the  signal  has  a  strong  periodic  component  with  a  period 
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FREQUENCY  IN  CYCLES  PER  SECOND 

Fig,-  12,  Power  spectral  density  estimate  for  the  17.9  me 

signal  power  level  fluctuations  from  1421  to  1445 
PST  on  August  4,  I960,  The  solid  and  broken  lines 
represent  two  separate  estimates  of  the  same  data. 
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of  a  day.  However,  this  ie  not  of  prime  interest  in  this 
study. 

The  power  spectral  density  estimates  for  the  shorter 
records  are  more  irregular  in  shape  than  those  for  the  56 
minute  transmission  period  of  WWV*  This  is  due  to  the 
length  of  the  data  samples  involved,  the  shorter  data 
samples  yielding  the  less  stable  estimates. 

Two  of  the  graphs,  Figures  11  and  12,  show  two  curves 
each,  one  solid-line  curve  and  one  broken-line  curve.  These 
two  curves  are  two  separate  estimates  for  the  same  data.  In 
the  case  of  Figure  11,  the  two  separate  estimates  were  made 
in  order  to  extend  the  range  of  the  estimate.  In  the  case 
of  Figure  12,  they  were  made  to  obtain  a  measure  of  the  ran¬ 
dom  errors  introduced  by  the  computing  equipment.  This  error 
lies  in  the  neighborhood  of  three  to  ten  percent  as  can  be 
seen  by  comparing  the  two  curves. 

The  estimate,  Figure  13,  for  the  day  of  September  22, 
i960,  from  1644  to  1658  PST  is  likely  to  be  of  low  validity. 
The  recorded  signal  was  of  such  small  amplitude  that  the  re¬ 
corder  noise  level  is  fairly  large  with  respect  to  the  sig¬ 
nal  level. 

Figures  11  and  12  show  two  estimates  for  data  samples 
that  are  almost  consecutive  in  time.  These  two  estimates 


peak  at  about  the  same  point  and  have  the  same  general 
shape  except  that  Figure  11  shows  a  greater  increase  in 
power  density  at  fairly  low  frequencies. 

Figures  l4  and  15  also  show  estimates  for  data 
which  are  consecutive  in  time.  These  two  estimates  drop 
off  to  zero  at  about  the  same  point  with  about  the  same 
rate.  Since  these  two  estimates  are  both  for  56  minute 
long  data  samples  they  should  be  reasonably  stable;  how¬ 
ever,  the  variance  of  the  estimates  is  not  known. 

One  interesting  facet  of  the  power  spectral  density 
estimates  is  that  the  frequency  range  of  the  fades  for 
the  fshorter  transmission  path  which  entails  mostly  one 
or  two  hops  (Figures  14  and  15)  is  almost  identical  with 
the  frequency  range  of  the  fadeB  for  the  longer  trans¬ 
mission  path  which  entails  one,  two,  threo  and  possibly 
four  hops  (Figures  10  through  15) •  If  the  fades  are  due 
mostly  to  beating  between  Doppler-shifted  signals  arriving 
via  different  modes  of  transmission  one  might  expect  a 
possibly  higher  frequency  range  for  the  greater  number  of 
possible  hops  since  the  signal  will  be  shifted  in  frequency 
at  each  reflection  from  a  moving  ionized  layer.  These  fre¬ 
quency  shifts  should  be  cumulative  at  least  part  of  the 
time. 


Since  the  quantity  of  d^ta  analyzed  is  still  small  the 
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certainty  of  any  conclusion  drawn  from  the  identical  fre¬ 
quency  ranges  is  poor* 

11,  Conclusions 

In  general,  the  fading  process  is  a  random,  quaBi- 
periodic  process  with  a  power  spectral  density  falling 
below  one  cps.  The  upper  limit  of  this  power  density 
seems  to  be  remarkably  constant  in  frequency. 

The  most  productive  method  of  estimating  the  power 
spectral  density  is  with  the  use  of  the  analog  computer. 
The  expected  value  of  the  analog  computer  estimates  will 
converge  to  the  true  power  spectral  density  for  higher 
and  higher  Q  of  the  simulated  resonant  circuit.  However, 
the  possible  spread  of  the  estimates  is  unknown.  It 
would  be  desirable  to  compute  the  possible  spread  of  the 
estimate  in  order  to  obtain  a  statistical  measure  of  cer¬ 
tainty  for  the  estimates.  The  mathematical  complexity 
of  this  computation  precluded  its  inclusion  in  this  re¬ 


port  . 
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APPENDIX  A  -  Equipment 

The  records  of  the  variation  of  the  amplitude  of 
the  incoming  wave  were  in  the  form  of  magnetic  tape  re¬ 
cordings  of  the  diode  detector  voltage  out  of  a  receiver 
type  P-390A/URR,  a  Collins  Radio  Company  receiver  type 
51 J4,  serial  number  784,  and  a  Collins  Radio  Company 
receiver  type  75^1,  serial  number  1524.  The  rf  gain  of 
the  receivers  was  adjusted  so  that  the  diode  output  was 
proportional  to  the  input. 

The  magnetic  tape  recorder  used  in  the  experiment 
was  a  Precision  Instruments  model  200  recorder  and  the 
tape  recordings  were  made  of  an  FM  signal  whose  frequency 
deviation  was  proportional  to  the  amplitude  of  the  sig¬ 
nal  within  plus  or  minus  one  percent*  The  primary  data 
are  stored  on  magnetic  tapes. 

The  tape  speedB  used  for  the  recording  were  7-1/2 
inches  per  second  and  0.075  inches  per  second  and  the 
harmonic  distortion  of  the  recording  system  is  less  than 
1-1/2%  from  0  to  1250  cps.  The  original  data  are  in  the 
form  of  records  lasting  from  fifteen  minutes  up  to  an 
hour.  For  purposes  of  computing  with  these  data  on  the 
digital  computer  the  original  tape  was  played  back  through 
a  direct  record  system  and  the  frequency  of  the  signal  was 


counted  with  a  Hewlett  Packard  counter  type  522B,  and 
then  printed  out  on  paper  as  a  sequence  of  numbers  with 
a  Hewlett  Packard  type  56OA  digital  recorder*  In  this 
fashion  the  signal  was  sampled  at  a  maximum  rate  of  one 
sample  every  one-tenth  of  a  second*  In  cases  where  the 
variation  of  the  data  was  small  over  periods  of  time 
that  were  large  with  respect  to  the  sample  interval  time, 
the  time  between  samples  was  increased  to  1/2  of  a  second. 

Directly  counting  the  frequency  of  the  PM  signal  on 
the  tape  introduces  an  offset  and  a  scale  factor  into 
the  original  signal  by  virtue  of  the  fact  that  aero  diode 
voltage  corresponds  to  6750  cps.  Plus  or  minus  b0%  fre¬ 
quency  deviation  from  this  depends  on  d.c.  amplifier 
gain.  Consequently  the  digital  computer  had  to  be  in¬ 
structed  to  convert  this  set  of  numbers  back  to  the 
original  data. 

Data  were  also  analyzed  by  using  the  analog  computer 
as  a  spectrograph.  To  speed  up  this  operation  yb  hod 
change  the  time  scale  of  the  records.  This  was  accom¬ 
plished  by  using  two  tape  recorders  together  and  by  re¬ 
producing  the  signal  at  a  high  tape  speed  on  one  machine 
and  re-recording  it  at  a  low  tape  speed  on  the  other 
machine.  The  total  time  reduction  obtained  by  the  pro¬ 
cess  was  52  or  64  times  depending  on  the  particular  set 
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of  data. 

The  second  tape  machine  UBed  in  this  process  was  an 
Ampex  model  tape  recorder.  The  shortened  recordings  ob¬ 
tained  with  this  process  were  then  used  in  con junction 
with  a  Berkeley  Ease  analog  computer. 

The  antenna  used  with  the  receiver  was  a  three  ele¬ 
ment  Yagi  antenna  mounted  approximately  sixty  feet  above 
the  ground.  One  of  the  transmission  distances  involved 
in  the  experiment  was  6000  km,  the  transmitter  being 
located  in  a  direction  approximately  14  degrees  south  of 
east  from  Seattle. 

The  transmitter  used  in  part  of  the  experiment  was 
transmitting  a  cw  signal  with  no  modulation  except  on 
the  hour  and  half-hour  times,  during  which  it  was  keyed 
on  and  off  for  identification  purposes.  The  other  trans¬ 
mitter  was  WYiTV  in  Belleville,  Maryland, 
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